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ABSTRACT 

A local luminosity function at 15/im is derived using the bivariate (15/xm vs. 60/xm 
luminosity) method, based on the newly published ISOCAM LW3-band (15/xm) survey 
of the very deep IRAS 60/xm sample in the north ecliptic pole region (NEPR). New 
IRAS 60/tm fluxes are obtained using the SCANPI/SUPERSCANPI software at the 
new ISOCAM positions of the sources in the sample. It is found to be in excellent 
agreement with the 15/im local luminosity function published by Xu et al (1998), 
which is derived from the predicted 15/tm luminosities of a sample of IRAS 25/tm 
selected galaxies. Model predictions of number counts and redshift distributions based 
on the local luminosity function and assumptions of its evolution with the redshift are 
calculated and compared with the data of ISOCAM 15/im surveys. Strong luminosity 
evolution on the order of L oc (1 + z) 4 ' 5 is suggested in these comparisons, while pure 
density evolution can be ruled out with high confidence. The sharp peak at about 
0.4mJy in the Euclidean normalized differential counts at 15/zm can be explained by the 
effects of MIR broadband emission features, eliminating the need for any hypothesis for 
a 'new population'. It is found that the contribution from the population represented 
by ISOCAM 15/im sources can account for the entire IR/submm background, leaving 
little room for any missing 'new population' which can be significant energy sources of 
the IR/submm sky. 



Subject headings: galaxies: luminosity function, mass function - galaxies: photometry 
- galaxies: starburst - galaxies: statistics - infrared: galaxies 



1. Introduction 

Our understanding of galaxy formation and evolution in the early epochs of the Universe has 
been vastly improved in the past a few years, thanks mainly to new deep surveys in a wide range 
of wavebands, ranging from the HST's WFPC2 (UV and optical) and NICMOS (NIR) surveys 
in the Northern and Southern Hubble Deep Fields (Williams et al. 1996; Williams et al. 1998; 
Thompson et al. 1999) to the SCUBA submm surveys (see Blain et al. 1999b for a review). In 
particular, several mid-infrared (MIR) and far-infrared (FIR) deep surveys (see Elbaz et al. 1998b, 
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Puget & Lagache 1998 for reviews) were conducted using the Infrared Space Observatory (ISO) 
(Kessler et al. 1996) during its 29 months mission (Nov. 1995 - April. 1998). The new results 
from these surveys (Aussel et al. 1999; Puget et al. 1999; Dole et al. 1999) indicate significant 
improvements in sensitivity and accuracy over the earlier published results (Rowan-Robinson et al. 
1997; Kawara et al. 1998). Strong cosmic evolution in the population of infrared-emitting galaxies 
is indicated in these results (Rowan-Robinson et al. 1997; Kawara et al. 1998; Aussel et al. 1999; 
Puget et al. 1999), consistent with the results of SCUBA surveys (Blain et al. 1999a) and with 
the scenario hinted at by the newly discovered cosmic IR background (CIB) (Puget et al. 1996; 
Hauser et al. 1998; Dwek et al. 1998; Fixsen et al. 1998). These results challenge those from 
the UV/optical surveys (Madau et al. 1998; Pozzetti et al. 1998) in the sense that substantially 
more (i.e. a factor of 3 - 5) star formation in the earlier Universe with z ^ 2 may be hinted at in 
the IR/submm counts and in the CIB (see, e.g. Rowan-Robinson et al. 1997) compared to that 
derived from the UV/optical surveys (Madau et al. 1998; Pozzetti et al. 1998). The reason of 
this discrepancy is attributed to dust extinction which may hide much of the star formation in the 
early Universe from the UV/optical surveys (see Lonsdale 2000 for a review). 

The best observed band in these ISO surveys is the ISOCAM LW3 (15/um) band, with 14 
surveys covering a wide range of flux density from 0.05 mJy to 50 mJy (Elbaz et al. 1998b). 
Compared to the longer wavelength ISOPHOT surveys (e.g. 175/xm FIRBACK survey, Puget 
et al. 1999), ISOCAM LW3 surveys have the advantage of using large detector arrays (32 x 32 
compared to the 2x2 array of the ISOPHOT-C200 camera) and having much better angular 
resolution which allowed the surveys to go very deep (~ 0.1 mJy compared to sensitivity limits of 
175/Um surveys of ~ 100 mJy) before reaching the confusion limit. In addition to the indication 
of the significant cosmic evolution of IR galaxies, two other very interesting results are emerging 
from these MIR surveys: 

1. The Euclidean normalized differential number counts of the 15/wn band have a sharp peak at 
about 0.4 mJy (Elbaz et al. 1999) which, Elbaz et al. (1998a) claim, can only be explained 
by adding to the number count model a new population of objects which emerge (with 
increasing z) rapidly after z ~ 0.4 and start to dominate the counts below 1 mJy, but 
contribute negligibly in brighter flux range (Fig. 10 of Elbaz et al. 1998a). 

2. Subject to the substantial uncertainties in their IR SEDs, the integrated light of the galaxies 
detected in 15 fj,m surveys may account for most the IR/submm background (Elbaz et al. 
1999). 

Taken at face values these results may have far-reaching impact on the studies of galaxy evolution, 
suggesting that the objects most responsible for the CIB, which are mostly missed by the 
UV/optical surveys, are already identified in the 15^m surveys and, perhaps more interestingly, 
they are from a new population not seen in the local Universe (i.e. not among the IRAS sources 
which in general have z ;$ 0.1). 
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However, scrutiny of these interpretations of the 15/xm counts is imperative because of the 
following two complications: (1) the effects of prominent emission features in the wavelength 
range of 3 - 13/xm (Puget and Leger 1989) which can cause very significant K-corrections in MIR 
surveys (Xu et al. 1998); and (2) the lack of a local luminosity function (LLF) in the 15/j.m band 
which is needed for the quantitative determination of the evolution rate from the number counts 
(to date, the IRAS \2fxm luminosity functions have been mostly used in the interpretation of the 
ISOCAM 15/xm counts, resulting in large uncertainties due to the significant variations of the 
fibaml '/i2«m among galaxies; see Elbaz et al. 1999). Indeed, the source count model of Xu et al. 
(1998), which takes into account the effect of the MIR emission features, did predict the bumps 
and dips in the counts similar to what is seen in the 15/xm surveys. Accordingly, Xu et al. (1998) 
made the warning that determinations of the evolution rate based on the slope of source counts 
will have to treat the effects of the MIR emission features carefully. Xu et al. (1998) also derived 
a LLF at 15/im from the predicted 15/xm flux densities of a sample of 1406 IRAS 25/xm selected 
galaxies based on a three component (cirrus/PDR, starburst, AGN) SED model. 

Recently, Aussel et al. (2000) published an ISOCAM 15/im survey of the very deep IRAS 
QOfim sample in the northern ecliptic pole region (NEPR) (Hacking and Houck 1987, hereafter 
HH87; Hacking 1987). In this paper, we will derive a new 15/im LLF based on the data of this 
survey, and on the new IRAS SCANPI measurements at 60/mt of the same sources using the new 
ISOCAM positions. This is then compared to the LLF of Xu et al. (1998). With the confidence 
in the 15/j.m LLF gained from this new study we explore further, using the model of Xu et al. 
(1998), the quantitative interpretation of results of the ISOCAM deep surveys as published by 
Elbaz et al. (1999) and by Aussel et al. (1999). 

2. New Ihum LLF from NEPR Sample 

2.1. New IRAS SCANPI measurements at ISOCAM positions 

The ISOCAM observations at 15/jm for 94 out of 98 galaxies in the very deep IRAS 60/j.m 
sample in the NEPR (HH87) are described out by Aussel et al. (2000). Altogether 106 sources 
were detected with signal to noise ratios > 3. Several IRAS sources correspond to multiple 
ISOCAM sources, given the much better angular resolution of ISOCAM (~ 10") compared to 
the IRAS resolution (> 1'). The mean position offset of the ISOCAM sources relative to IRAS 
positions (HH87) is ~ 10". 

New IRAS measurements at 60/um using the IPAC software SCANPI were carried out at 
the positions of the 117 (11 having signal to noise ratio < 3) ISOCAM sources in the NEPR 
sample listed in the Table 3 of Aussel et al. (2000), exploiting the IRAS survey database. This 
is to get a better corresponding 60fim flux density for each ISOCAM source, and also to bring 
the IRAS flux densities obtained by HH87 to the newest IRAS standard (Moshir et al. 1992). 
For this task, three SCANPI queries were run and finished on Oct. 28, 1999, using the default 
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Fig. 1. — Comparisons between feoum from this work (new SCANPI/SUPERSCANPI reductions) 
and from HH87 (Fig. la), and from FSS (Fig. lb). 
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settings of SCANPI 1 . Of the 117 sources, 99 are detected at 60/im. For the 18 ISOCAM sources 
undetected by SCANPI, the interactive SUPERSCANPI has been run in attempts to increase the 
sensitivity when, in addition to the survey data, data from pointed observations (HH87; Gregorich 
et al. 1995) are also available. As in HH87 we include only pointed observations with 'deep-sky' 
macros in the SUPERSCANPI processing. Five additional detections were obtained through the 
SUPERSCANPI coadds. This leaves 13 ISOCAM sources undetected in the 60^m band, which 
have in general large offsets (<offset>~ 60") from the original IRAS positions. Of the 104 sources 
detected by SCANPI and SUPERSCANPI in the 60/xm band, many (usually those corresponding 
to the same source in the HH87 list) are confused. The plots of SCANPI and SUPERSCANPI 
processing of these confused sources awere manually inspected to determine the best total flux 
densities of these confused sources, which are assigned to them jointly. A final list of 106 sources 
is given in Table 1, including the sources undetected by SCANPI/SUPERSCANPI, and the 4 
sources for which Aussel et al. reported only upper-limits at 15/xm. Redshifts are found for 68 of 
them from Ashby et al. (1996). 

The SCANPI/SUPERSCANPI results are compared with the flux densities reported by 
HH87 and with those from the IRAS Faint Source Survey (FSS) (Moshir et al. 1992) for sources 
which observe the following criteria: (1) ISOCAM position within 20" of the IRAS position, (2) 
not confused with any other sources in the IRAS 60/im band, and (3) for FSS comparison, they 
have to be listed in the Faint Source Catelog (FSC) or in the Faint Source Reject File (FSR) in 
the IRAS database. Thirty sources are selected by the first two criteria, 27 of them also pass 
the criterion (3). It is found (Fig. la and lb) that the 60 \im flux densities obtained in this paper 
are consistent with those listed in the FSC and FSR, but about 20% higher than those reported 
by HH87. Given that the FSC and FSR results represent the new standard of IRAS products, 
it is likely that the lower flux densities of HH87 are due to some systematic biases in the early 
processing of IRAS data. 

2.2. New 15/xm LLF Derived from Bivariate 15/xm/60/xm LF 

Exploiting the NEPR sample, a LLF at 15/xm can be constructed using the so called 'bivariate 
method' (see, e.g., Xu et al. 1998), transferring the 60/wn LLF of IRAS galaxies, which has been 
well studied in the literature (Soifer et al. 1987; Saunders et al. 1991; Yahil et al. 1992), to 
15fim LLF utilizing the Li^ um / Leoum ratio v.s. L§§ U m relation. We include only the 64 sources in 
Table 1 which are detected in the 60/um band (including sources with upper-limits at 15/^m) and 
which have measured redshifts (Ashby et al. 1996). 

There are several concerns with regard to the sample: 



^or details of SCANPI and SUPERSCANPI processing, see the webpage 

http : / /www Apac.caltech.edu/ipac/iras/ scanpi-over.html. 
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Fig. 2. — NEPR sources (64) used in the bivariate analysis: Li5 Mm / L^o^m ratio versus 15/xm 
luminosity. 
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• The sample is incomplete. In particular, the requirement of redshift availability excludes 
about one third of the sources in Table 1. 

• Possible misidentification between the sources in the redshift survey of Ashby et al. (1996) 
and the sources in this work. 

• Given the depth, the sample is not really local (many sources have z > 0.1), hence may be 
affected by galaxy evolution with increasing redshift. 

• The redshift distribution of the sources shows strong clustering (Ashby et al. 1996). 

Will these affect the bivariate 15 [imj '60 fim luminosity function? The answer to this question 
depends on whether the 15iM7i-to-60//m color ratio is a sensitive function of the luminosity. This is 
because all the above potential problems with the sample are related to the redshift, and result in 
uncertainties in the luminosity distribution (the 'visibility function') but not in the £i5^m/-^60/wn 
ratio distribution. If the color ratio is insensitive to the luminosity, the conditional probability 
function 0(^15^/^60^1^60^™) (cf. Eq(7) of Xu et al. 1998), which converts the 60/jm LLF 
to the 15/xm LLF, will be rather constant and won't be affected significantly by uncertainties 
associated with the luminosity. 

Indeed, as plotted in Fig.2, the I'i5 //m /-^60/im ratio of sources in the NEPR sample appears 
to be rather insensitive to the luminosity. This result is similar to that of Soifer and Neugebauer 
(1991) who found that the -^25^m/-^60/im ratio does not depend on the infrared luminosity. 
This may not be surprising given that the mechanisms of the 15/j.m emission and of the 25/xm 
emission are nearly the same, namely due dominantly to the emission of small grains undergoing 
temperature fluctuations in normal galaxies such as the Milky Way, and to the warm dust emission 
associated with star formation regions in starburst galaxies such as M82 (Desert et al. 1990). This 
dualism in the radiation mechanism of the MIR continuum, at the wavelengths not contaminated 
by the MIR emission features, is the major reason for the lack of dependence of the two color ratios 
(L^m/Leo^m and L 2 5^ m / L mf , m ) on luminosity. At the same time, the very cold L 2 5^ m /Lm^m 
ratios of ultraluminouse galaxies (ULIRGs), which are due mostly to extinction at 25 fim (Xu 
and De Zotti 1989), also further weaken any statistical dependence of the L25 fim /L^o fim ratio on 
the grain temperature, the latter is a strong function of the luminosity as demonstrated by the 
Leo^m/Lioo/im v.s. luminosity relation (Soifer and Neugebauer 1991). 

The algorithm and the formulation used in this work are the same as presented in Section 3 
of Xu et al. (1998) 2 . The 60/xm LLF derived by Saunders et al. (1991) using the so-called 'non- 
parametric maximum-likelihood' method is taken, for which the effects of spatial galaxy density 
fluctuations, in particular the local over-density due to the local super cluster, are minimized 
(Saunders et al. 1991). In Table 2 the derived 15/xm LLF is listed, with Li^^ m being defined by 



2 There was an error in Eq(14) of Xu et al. (1998), which should have been Covar(Fi-i, Fi) — F t x Var(Fi-i)/ Fi-i. 
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Fig. 3. — 15fim local differential luminosity function. 
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vL v at 15pm and bin width <51og(Li5 Mm ) = 0.4. In Fig.3 this new 15pm LLF is compared to the 
15pm LLF of Xu et al. (1998) which is derived from the predicted 15pm luminosities of a 25pm 
selected sample of IRAS galaxies. Excellent agreement, in particular near the knee of the LLF 
(~ 10 95 L©), is found between these two LLFs, which are derived from completely different data 
sets using very different approachs. This verifies that both LLFs are reliable within the limits of 
their uncertainties. On the other hand, the two LLFs are complementary to each other. While the 
new LLF is derived from real 15 pm ISOCAM data obtained by Aussel et al. (2000), the size of 
this data set (64 galaxies) is much smaller than the IRAS sample (1406 galaxies) used by Xu et 
al. (1998). Consequently, the new LLF does not extend beyond L* as far as the LLF of Xu et al. 
(1998), namely being truncated at 10 113 L and with the point at 10 10,9 missing since there is no 
galaxy in that bin. It should be noted that both the 60pm LLF of Saunders et al. (1991) on which 
this work is based, and the 25pm LLF of Shupe et al. (1998) on which the 15pm LLF of Xu et 
al. (1998) is based, are derived using the maximum-likelihood method which minimizes the effect 
of density fluctuations. Also both the normalizations of the 60pm LLF of Saunders et al. (1991) 
and the 25pm LLF of Shupe et al. (1998) are carefully determined, and are transferred by the 
bivariate analyses to the 15pm LLFs of Xu et al. (1998) and of this work, respectively. The very 
good agreement between the points of the two 15pm LLFs near the knee (Fig.3) demonstrates 
that both normalizations are indeed reliable, eliminating a large uncertainty in the prediction of 
the local 15pm counts (Elbaz et al. 1999). In what follows we will use the LLF of Xu et al. (1998) 
when a 15pm LLF is needed. 



3. Galaxy Evolution Indicated by ISOCAM 15pm Counts 

In Fig. 4 we reproduce the results of ISOCAM 15pm surveys presented in Elbaz et al. (1999). 
In addition, counts derived based on the predicted 15pm flux densities of sources in the IRAS 
25pm selected sample of Xu et al. (1998) are plotted at flux density levels > 0.2 Jy. Counts 
fainter than 0.2 Jy are not plotted since they drop dramatically, indicating the increasingly severe 
incompleteness, due to the fact that the sample is 25pm selected rather than 15pm selected. Note 
that at the bright end (> 0.5Jy) the normalized counts of these sources are significantly higher 
than the counts in the fainter flux density bins. This excess of counts is very likely due to the 
overdensity associated with the local supercluster (Lonsdale et al. 1990; Saunders et al. 1991). 

As pointed out by Elbaz et al. (1999), the Euclidean normalized ISOCAM 15pm counts 
have a narrow and prominent peak at about 0.4 mJy. There have been suggestions that this 
peak indicates a new population of infrared sources emerging after redshift z ~ 0.4 (Elbaz et al. 
1998a). On the other hand Xu et al. (1998) argued that such a feature can be caused by the 
broad-band emission features often associated to the polycyclic aromatic hydrocarbon molecules 
(PAH features, see Puget and Leger 1989). However the model predictions published by Xu et al. 
(1998), specified by two galaxy evolution models including one pure luminosity evolution model 
with evolution rate as L = Lq x (1 + z) 3 and one pure density evolution with p = p$ x (1 + z) A , 
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Fig. 4. — Euclidean normalized 15/xm differential counts: model predictions compared to the 
observations. The observational data are taken from Elbaz et al. (1999). Data points (same as in 
Elbaz et al. 1999): A2390 (stars); ISOHDF-North (open circles), ISOHDF-South (filled circles), 
Marano FIRBACK Ultra-Deep (open squares), Marano Ultra-Deep (exes), Marano FIRBACK 
Deep (stars), Lockman Deep (open triangles), Lockman Shallow (filled triangles). At the bright 
end (> 0.2Jy) plotted (filled squares) are the counts derived based on the predicted lbfxm flux 
densities of the sources in the sample of Xu et al. (1998). 
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calculated using the number count model which takes into account the effects of these emission 
features, underestimated the counts when compared with the ISOCAM data (Elbaz et al. 1998a). 
This suggests that the evolution endured by the ISOCAM sources is stronger than that assumed 
by Xu et al. (1998) which is based on previous studies of IRAS counts (Lonsdale et al. 1990; 
Saunders et al. 1991; Pearson and Rowan- Robinson 1996). 

Here we present new model predictions using the same number count model of Xu et al. 
(1998), but with stronger evolution rates, and compare them with the ISOCAM 15^m data. We 
have assumed that galaxy formation starts at z=5 (the counts are not sensitive to this parameter). 
The cosmology adopted in the models plotted in Fig.4 is specified by Q m = 0.3 and = 0.7, 
which is suggested by recent observations of type I supernovae in distant galaxies (Garnavich et 
al. 1998). Models with £l m = 1 17 a = and with Q m = 0.3 Q\ = were also calculated, but 
not plotted here. The results from the f2 m = 0.3, S7a = cosmology are very close to the results 
presented in Fig.4, while the results from the Q m = 1, J7a = model fit the data slightly less well. 

The solid line represents the pure luminosity evolution model with Lex (1 + z) 4 - 5 and with 
a turnover at z = 1.5 beyond which the evolution turns flat (i.e. L=constant for z > 1.5). The 
short-dashed line gives the counts predicted by another pure luminosity evolution model with 
Lx (1 + z) 4 - 8 and with a turnover at smaller redshift: z=l (i.e. L=constant for z > 1). Finally 
the long-dashed line is the prediction for counts by a pure luminosity evolution model without 
any turnover, and with an evolution rate of Loc (1 + z) 4 ' 3 . Among the three luminosity evolution 
models, the one with turnover at z=1.5 (the solid curve) gives the best fit, closely reproducing 
the overall shape and the level of the observed counts. The model with a turnover redshift z=l 
predicts a peak which is too flat compared to the data. This is because, in the framework of the 
Xu et al. (1998) model, the broadband MIR features in the wavelength range 6 - 8.5//m, which 
are redshifted into the ISOCAM LW3 bandpass when z=l±0.2, are largely responsible for the 
narrow peak of the counts in Fig.4. When the turnover occurs at z=l, a significant number of 
sources in the redshift range of z=l - 1.2 are dropped compared to the models without turnover 
or with the turnover at 1.5, resulting in a less prominent peak. The model without turnover (the 
long-dashed line) predicts a peak at fi5 Mm ~ 0.2 mJy instead of 0.4 mJy as shown by the data. 

The dotted line shows the prediction by a density evolution model with comoving density 
p oc (1 + z) 7 - 5 until z=1.5, turning flat afterwards (i.e. /?=constant when z > 1.5). This model 
gives a reasonable fit to data points brighter than fi5 Mm ~ 0.2 mJy. However in the fainter flux 
levels, instead of turning down, the model prediction keeps rising in the plot until fi5^ m reaches 
as low as ~ 0.04 mJy. This is very different from the trend shown by the data. The reason for 
the difference between the density evolution model and the luminosity evolution models is that 
for a given redshift, say z=l which allows the 6 - 8.5/um emission features to be included in the 
LW3 bandpass, galaxies are much fainter in pure density evolution models compared to those in 
pure luminosity evolution models, therefore the bump caused by the K-correction due to the MIR 
emission features occurs in much fainter flux density levels (see Xu et al. 1998 for a more detailed 
discussion). 
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Fig. 5.— Redshift distribution of ISOCAM lh^m sources in the field of ISOHDF-North (Aussel et 
al. 1999): model predictions compared to the observational data. 
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Strong constraints on galaxy evolution can be obtained when the redshift information of a 
flux-limited sample is available. Aussel et al. (1999) found in the literature 29 redshifts for 49 
sources in the main source list of the ISOCAM HDF (North) survey. In Fig. 5 the histogram of the 
redshifts of 17 of these sources with fis^m > 0.1 mJy is compared to the model predictions by the 
three models presented in Fig.4. A sky coverage of 16 acrmin 2 (Aussel et al 1999) and a correction 
factor of 2 for the incompleteness are assumed in these calculations. All three luminosity evolution 
models give reasonably good fits to the data in the bins with z < 1, and over-predict the counts 
in bins of z > 1. In particular, both the model with turnover at z=1.5 and the model without 
turnover predict some sources (6 by the former and 10 by the latter) with z> 1.4 while none are 
found in the data. For the model with turnover at z=1.5, the best fitting model in Fig.4, the 
missing of sources at z > 1.4 could be due to small number statistics or to the incompleteness 
of the data at high redshifts. This highlights the demand for larger and more complete redshift 
samples of ISO sources. In fact, when redshifts for a large (a few hundred) and complete flux 
limited sample are available, luminosity functions of ISO sources can be calculated for different 
redshift epochs which will give the most direct information about the evolution of these sources. 

The reasonably good fits to both the source counts (Fig.4) and the redshift distribution 
(Fig. 5) by the three luminosity evolution models (in particular the model of Loc (1 + z) 4 ' 5 with a 
turnover at z = 1.5) demonstrate that indeed the narrow peak of the ISOCAM 15/xm counts can 
be well explained by the effect of broadband MIR emission features which is the essential element 
of the model of Xu et al. (1998), and there is no need to invoke a 'new population'. At the same 
time, the evolution rate implied by the model fit is much stronger than those given by previous 
studies on IRAS sources (Loc (1 + z) 3 , Pearson and Rowan-Robinson 1996; Lonsdale et al. 1990), 
but is consistent with what is found in the UV and optical deep surveys (Lex (1 + z) s - 95±0 - 75 , Lilly 
et al. 1996). Analyzing multiband data from IRAS, ISO, SCUBA and COBE, Blain et al. (1999a) 
also obtained a relatively high evolution rate ((Loc (1 + z) 3 - 8±0 2 ). Given the lack of dependence of 
the Lis^m/Lso^m ratio on the luminosity, one expects that the evolution rates of the L\^^ m and of 
the 

L*60fj, m should be similar. In Fig. 6 we compare the counts at QOfim predicted by a luminosity 
evolution model assuming L,Qo^ m oc (1 + z) which turns flat {Lqq^„ l =constant) at z = 1.5, with 
IRAS data. The large filled circles are the counts from this work (new NEPR sample), which 
are about 30% (i.e. ~ 0.12 dex) higher than those of HH87 (crosses). This discrepancy can be 
fully explained by the fact that the 60/xm fluxes obtained by the new SCANPI/SUPERSCANPI 
processing are about 20% higher than those of HH87 (Fig. la), given that the Euclidean normalized 
differential counts scale with the flux to the 1.5 power. At the same time, the model predictions 
(solid line) indeed reproduce the trend and the level of those more recent IRAS counts (Bertin et 
al. 1997; Gregorich et al. 1995; Saunders et al. 1991; Lonsdale et al. 1990; Rowan-Robinson et al. 
1990) quite well. 

Pure density evolution, which gives a poor fit to the number counts, can be ruled out with 
high confidence. Given that any density evolution will push the peak in the number counts to 
fainter flux levels than shown by the data, it is quite certain that, as far as the choice between 
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Fig. 6. — Euclidean normalized differential counts of the IRAS 60fim band: model predictions 
compared to the observations. Data points: new NEPR (this work) (large filled circles); HH87 
(exes); Gregorich et al. 1995 (open stars); Bertin et al. 1997 (open circles); Lonsdale et al. 1990 
(filled squares); Saunders et al. 1991 (filled squares), Rowan- Robinson et al. 1990 (open squares). 
Model: LQQ^ m oc (1 + z) 4 ' 5 when z < 1.5, -Z^o^m =constant when z > 1.5 (solid line). 
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density and luminosity evolution is concerned, luminosity evolution is the dominant cause of the 
high number counts. This conclusion is in agreement with that of Blain et al. (1999a) which is 
obtained from a completely different argument, namely that a pure density evolution model which 
can fit the IR/submm counts will produce too much IR/submm background. 

Our results show that the narrow peak of the ISOCAM 15/um counts at about 0.4 mJy may 
not be used as an evidence for a 'new population' of faint MIR sources. On the other hand, 
a luminosity evolution in the luminosity function of infrared galaxies, as suggested by our best 
fitting model, does not necessarily mean that it is the same galaxies that we are seeing in the local 
Universe that are shining tens or even hundreds times brighter in the early epochs of the Universe. 
Indeed, the preliminary results of optical identifications of ISOCAM LW3 sources indicate that 
beyond z ~ 0.7 most of them are interacting galaxies (Aussel et al. 1999; Elbaz et al. 1999), 
while the local MIR selected extragalactic sources are mostly single late type galaxies similar to 
the Milky Way (Rush et al. 1993). Given the high incompleteness of the redshift data, at this 
stage the major constraint on the evolution of ISOCAM sources is from the counts, which is 
mostly determined by how the luminosity function evolves around L* (the luminosity distribution 
of sources in a given fis^m bin peaks strongly around L*). Therefore, without any extrapolation, 
what we know now is that whatever the population of IR sources is at redshift z ~ 1, its comoving 
density is about the same as the IR galaxies in the local Universe, while the characteristic 15fim 
luminosity of the faint sub-mJy sources is about 20 times the L* of local IR galaxies, namely at 
the ~ 10 11 L level. If these sources are indeed similar to the local gas-rich spiral-spiral galaxy 
pair systems, which dominate the bright end (Lfi r > 2 x 10 11 L Q ) of IRAS luminosity function (Xu 
and Sulentic 1991), the implied density enhancement of these sources at z ~ 1 compared to their 
density in the local Universe is more than an order of magnitude. Although the population of 
these interacting galaxies is not really 'new' (i.e. they are already important contributors of MIR 
counts in the local Universe), it is quite possible that this population may evolve much faster than 
normal late-type galaxies, and even than AGNs. For the sake of simplicity, we have treated all IR 
galaxies as a single population in our model and have not considered any 'differential evolution' 
(i.e. different evolution rates for galaxies with different luminosities). When more constraints on 
the nature of 15fim sources are available from future follow-up observations, a model treating 
the evolution of different galaxy populations differently (e.g. separating interacting galaxies 
from single galaxies and AGNs), such as the model by Franceschini et al. (1994), will be more 
appropriate. 

Our results also suggest that the MIR emission features are present in the SEDs of galaxies 
with redshifts up to z ~ 1. Whether this is still true for galaxies with even larger redshifts will be 
found out by the future SIRTF mission (Cruikshank and Werner 1997). If so, these features will 
facilitate a powerful new method of obtaining redshifts in infrared for the optically faint, heavily 
extinguished galaxies. 

As for the question of whether galaxy evolution has a turnover at z=l-2, our results are not 
conclusive, though a positive answer is favoured by the model fits (Fig.4). Given the significant 
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effect of the MIR emission features, which happens to affect the 15/xm counts at redshift just 
below z=1.5, the turnover favoured by our model may well be a false signal. Deep surveys at 
longer wavelengths (e.g. at 25 fxm and 70fim using the SIRTF/MIPS detector arrays) where the 
MIR features will be redshifted into the bandpass at larger z, are certainly desirable for the 
determination of the evolution of IR galaxies beyond z=1.5 (Xu et al. 1998). 

How much do these sources contribute to the Cosmic IR/sub-mm Background radiation? To 
answer this question, a Monte Carlo simulation based on the source count model of Xu et al. 

(1998) was carried out, assuming the galaxy evolution model that gives the best fit to the data in 
Fig.4 (Loc (1 + z) 4 ' 5 with a turnover at z = 1.5). Sources with certain Li^^ m and z are generated 
according to predictions of the number count model, and IR SEDs taken from the SED sample 
of 1406 galaxies in Xu et al. (1998) are assigned to these sources in accordance with their rest 
frame L\^ m . It should be noted that the SEDs modeled by Xu et al. (1998) stop at 120^m. In 
order to estimate the contribution from MIR galaxies to the IR background radiation at longer 
wavelengths, we have assumed that the IR emission at A > 120^m of all sources generated in 
the simulation has the same spectrum specified by a modified blackbody with T = 40K and the 
emissivity index of (3 = 1.5 (Blain et al. 1999). In reality, this submm SED may only apply to 
the luminous IR starburst galaxies (LIRGs) while the SEDs of less active galaxies are likely to be 
much colder (Eales et al. 1999). However, since the largest contribution to the submm background 
is from LIRGs (Blain et al. 1999), we neglect this complication here. 

For any given wavelength, all the fluxes from these simulated sources are summed up, 
resulting in the predicted contribution of the population of 15/xm sources at the wavelength in 
question. Again we assume that galaxy formation starts at z=5. In Fig.7, this prediction (the 
solid curve) is plotted against several measurements/upper- limits of the IR/submm background. 
The upper-limits are all from the studies of TeV gamma-ray sources (Dwek &: Slavin 1994; Stanev 
& Franceschini 1998). The filled circles with error bars are COBE/DIRBE results taken from 
Lagache et al. (1998), and the large crosses are from SCUBA results (Blain et al. 1999). The two 
dashed curves outline the range of the submm IR background detected by COBE/FIRAS (Fixen 
et al. 1998). According to our results, the contribution of the MIR galaxies to the background 
between 10-30/um is at the 4 nW/m 2 /sr level, compared to the results reported by Elbaz et al. 

(1999) that the 15fj,m background due to sources brighter than 50^Jy is 3.3 nW/m 2 /sr. Note that 
this already meets the upper-limits obtained by Stanev and Franceschini (1998) from the analysis 
of TeV emission of Mrk501. At longer wavelengths, the predicted contribution to the background 
emission agrees very well with the DIRBE points, and lies slightly above the upper-boundary of 
the measured submm background. Compared to the results of previous calculations on the cosmic 
IR background using 'backward evolution' models (Hacking and Soifer 1991; Beichman and Helou 
1991; Malkan and Stecker 1998), the category this work belongs to, our result is about a factor of 
2 higher because the evolution rate hinted at by ISOCAM 15/t/m surveys is significantly stronger 
than those used in the previous works. Our result is in agreement with Elbaz et al (1999), who 
found that the ISOCAM 15/xm sources may be able to account for the majority of the IR/submm 
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Fig. 7. — Cosmic IR/submm background: model predictions compared to the observational data. 
Symbols: Model prediction (z < 1.5: Loc (1 + z) 4 ' 5 ; 1.5 < z < 5: L=constant) — solid curve; 
COBE/DIRBE results (Lagache et al. 1998) — large rilled circles with error bars; SCUBA results 
(Blain et al. 1999) — large crosses; the range of COBE/FIRAS results (Fixen et al. 1998) — two 
dashed curves; upperlimits from TeV gamma-ray radiation of Mrk403 and Mrk501 (Dwek & Slavin 
1994; Stanev & Franceschini 1998) — diamonds and exes with upper-limits. 
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background. Taken at face value, the result in Fig.7 indicates that nearly all of the sources 
contributing significantly to the IR/submm background are already present in the population of 
15/xm sources detected by ISO, and very little room is left for any missing 'new population' which 
can be significant energy sources of the IR/submm sky. 
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Table 1. New feo^m for the NEPR sample 
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Table 1 — Continued 
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86.9 





Table 1 — Continued 



LP. 1 f 60Atm 2 error 3 fgff m 4 offset 5 redshift 6 



(HH87) 


(mJy) 


(mJy) 


(mJy) 


(") 




3-89 


90 


16 


69 


56.5 


0.029 


3-90 


100 


9 


65 


8.6 


0.0720 


3-91 


80 


11 


58 


15.2 




3-92abc 


90 


19 


75 


56.9 


0.0880 


3-93 


140 


19 


99 


23.0 


0.0690 


3-94a 


130 


16 


98 


3.7 


0.0502 


3-94b 




20 




89.3 




3-96 


190 


18 


150 


5.8 


0.0321 


2-16 


40 


13 


38 


9.5 




3-06 






62 




0.0754 


3-22 






88 






3-87 






55 






3-95 






78 




0.187 



These are the same sources as listed in Table 3 of Aussel 
et al. (2000). Sources confused with each other in the 60/um 
band are grouped together as single entries. 

2 60/im flux density obtained in this work. For 
undetected sources the entry is blank. 

3 The one-<7 error of fm^m- 

4 60/im flux density taken from HH87. When more than 
one source corresponds to a single HH87 source, the flux 
density is assigned to the one with the least offset. 

5 Offset of the 15/xm source from the HH87 source, 
taken from Aussel et al. (2000). When more than one 
15/im source corresponds to a single source here (confused 
sources), the smallest offset is taken. 

6 Redshift taken from Ashby et al. (1996). 
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Table 2. l5/j,m Local Luminosity Function (LLF) from the NEPR sample 



\og{vL v {lhnm)/L & ) 


log(0/(Mpc 3 mag 1 )) 


1 a error 


8.1 


-3.70 


-3.71 


8.5 


-2.42 


-2.78 


8.9 


-2.68 


-3.40 


9.3 


-3.09 


-3.56 


9.7 


-3.47 


-4.21 


10.1 


-4.11 


-4.43 


10.5 


-5.08 


-5.64 


11.3 


-7.49 


-7.49 



